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Abstract. For a paradigmatic model of chemotaxis, we analyze the effect of how 
a nonzero affinity driving receptors out of equilibrium affects sensitivity. This 
affinity arises whenever changes in receptor activity involve ATP hydrolysis. The 
sensitivity integrated over a ligand concentration range is shown to be enhanced 
by the affinity, providing a measure of how much energy consumption improves 
sensing. With this integrated sensitivity we can establish an intriguing analogy 
between sensing with nonequilibrium receptors and kinetic proofreading: the increase 
in integrated sensitivity is equivalent to the decrease of the error in kinetic 
proofreading. The influence of the occupancy of the receptor on the phosphorylation 
and dephosphorylation reaction rates is shown to be crucial for the relation between 
integrated sensitivity and affinity. This influence can even lead to a regime where 
a nonzero affinity decreases the integrated sensitivity, which corresponds to anti¬ 
proofreading. 


1. Introduction 

Bacterial chemotaxis, a process by which the cell directs its motion in response to 
external ligand concentration, is a canonical example of biological sensing. Experiments 
with E. coli have provided much insight into chemotaxis [1,2], making this bacterium 
sensory system a particularly well understood example. In E. coli chemotaxis, the 
sensitivity is a key observable quantifying the response in activity inside the cell due to 
a change in the external ligand concentration. 

Stochastic models for E. coli [3-6] receptors often assume that changes in activity 
are described by an equilibrium process involving only conformational changes, leading 
to an equilibrium Monod-Wyman-Changeux (MWC) model [7,8]. However, chemical 
reactions where the receptor changes from an inactive to an active state often involve free 
energy consumption through, for example, adenosine triphosphate (ATP) hydrolysis. A 
stochastic model including this feature must have transition rates that break detailed 
balance leading to a nonzero affinity, corresponding to the chemical potential difference 
involved in ATP hydrolyisis, driving the process out of equilibrium. 
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Two recent studies have analyzed the effect of such an affinity in models related to 
the E. coli sensory network. Tu [9] has considered the effect of the driving affinity on 
both the dwell-time distribution and the sensitivity in a model for the flagellar motor 
switching between run and tumble. Skoge et ah [10] have shown that nonequilibrium 
receptors can increase the signal to noise ratio for fixed sensitivity. 

Beyond E. coli chemotaxis, the effect of energy dissipation in biological processes 
involving information processing has received much attention recently [11-21]. A 
prominent example among such processes is kinetic proofreading [22-25], which is a 
dissipative error reduction mechanism related to copying biochemical information. As 
this error reduction is achieved through free energy consumption, a nonzero affinity 
driving the process out of equilibrium is also present in kinetic proofreading. Specifically, 
relations between the error and the driving affinity have been obtained [25, 26] (see 
also [27-32] for other recent works). 

In this paper, we consider a nonequilibrium model for E. coli receptors including 
ATP hydrolysis in the chemical reactions that involve changes in activity. We quantify 
the effect of having a nonzero driving affinity on sensing by analyzing an integrated 
sensitivity, which is an integral of the sensitivity over a concentration range. This 
observable is shown to have a simple relation with the affinity driving the process 
out of equilibrium. We show that sensing with nonequilibrium receptors and kinetic 
proofreading can be viewed as equivalent problems, with the increase in the integrated 
sensitivity in nonequilibrium sensing being analogous to the error reduction in kinetic 
proofreading. 

The transition rates for changes in activity are assumed to depend on whether 
the receptor is occupied by a ligand or empty. We show that this dependency is quite 
important for the relation between sensitivity and the driving affinity. There is even a 
regime where energy dissipation leads to a decrease in the integrated sensitivity, which 
is equivalent to an anti-proofreading regime in kinetic proofreading [32]. 

The paper is organized as follows. Section 2 contains a simple stochastic model for a 
single nonequilibrium receptor. In section 3 we introduce the integrated sensitivity and 
obtain its relation with the affinity driving the process out of equilibrium. The analogy 
between nonequilibrium sensing and proofreading is established in section 4. In section 
5, with a more general model for a single receptor, we analyze how the influence of the 
occupancy of the receptor on the phosphorylation and dephosphorylation reaction rates 
affects the relation between integrated sensitivity and affinity. We conclude in section 6. 
Moreover, Appendix A contains a generalization of the single receptor model analyzed 
in the main text to an arbitrary number of binding sites. 

2. Nonequilibrium receptor model 

The single receptor model we analyze in this paper is defined as follows (see Fig. 1). 
There are two binary variables a and b characterizing the state of the receptor, with 
6 = 1 if the receptor is occupied by a ligand (bound) and 5 = 0 if the receptor is free 
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Figure 1. Four-state model for a single receptor. Vertical transitions correspond to a 
change in activity, while horizontal transitions correspond to a change in the occupancy 
of the receptor. The phosphorylation rates in (4) are chosen as and 

= 7 ae"^^. The dephosphorylation rates in (5) are chosen as a;+ = aya{Ki/Ko)e^^ 
and uj- = ay a.. 


(unbound), and a = 1 if the receptor is active and a = 0 if the receptor is inactive. 

In equilibrium the free energy of the four different states can be written as [8,20] 

F{a,b) = aAE — b\n. —, (1) 

where AE is the conformational energy difference between active and inactive for a free 
receptor (6 = 0), Ka is the dissociation constant that depends on the activity a, and c is 
the external ligand concentration. Setting Boltzmann’s constant and the temperature 
to k-^T = 1, the equilibrium stationary probability is Fa^ oc exp[—F(a, 6)]. Denoting 
the coarse-grained probability by Pa = J2b' Pa,b', we obtain 


Pi 


= e 


AE 


eq 



( 2 ) 


where Pq = 1 — Pi. The average activity 


(ci)c — ^ ^ ^Pg 

a 


( 3 ) 


is just (a)c = El- It is assumed that the dependence of the dissociation constant on the 
activity is such that Ki > Kq, i.e., the free energy barrier for binding a ligand to an 
inactive receptor is smaller. Hence, from Eq. (2) the average activity is a decreasing 
function of the concentration. This single receptor MWC model already contains the 
key feature of self-regulation. However, in order to have cooperativity, which is another 
important feature of the MWC model, we need more than one binding site [8]. The 
generalization of the model to an arbitrary number of binding sites is contained in 
Appendix A. For our present purposes it is more convenient to restrict to a single 
binding site in the main text. 

We now consider a nonequilibrium model that includes ATP hydrolysis. For 
simplicity we assume that when the receptor is unbound only phosphorylation takes 
place and when the receptor is bound only dephosphorylation occurs. A more general 
model with phosphorylation and dephosphorylation occurring for both bound and 
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unbound states, and the implications of this generalization are discussed in section 
5. The phosphorylation reaction is represented as 


ATP + 


^^Ladp + J], 


where and denote transition rates. The dephosphorylation reaction reads 


/. }l 


( 4 ) 


( 5 ) 


where u:\ and uif are transition rates. 

With the free energy (1), the generalized detailed balance relation [33] imposes the 
following constraints on the rates. First, we have 


, K+cal Ki 


( 6 ) 


where Ap = /xatp — /^adp — /^P; is the free energy dissipated in one ATP hydrolysis. 
Second, the transition rates from 6 = 0 to 6 = 1, denoted by Wq^, and from 6 = 1 to 
6 = 0, denoted fulhll the relation 


w, 


\n^ = \n{c/Ka). 


w 


( 7 ) 
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With these two constraints the product of the transition rates in a cycle in the clockwise 
direction in Fig. 1 is precisely A/r, which is the affinity driving the process out of 
equilibrium. For simplicity we use the specihc transition rates given in Fig. 1. The 
parameters 7a and 7b set the time-scale of the active/inactive and bound/unbound 
transitions, respectively. The parameter a is related to redistributing energy weights 
among the transition rates in such way that the constraints (6) and (7) are still fulhlled. 

A reasonable assumption is that ligand binding is much faster than activity changes, 
i.e., 7a/7b 1 [3]. With this assumption, calculating the stationary probability 

distribution using standard methods [34,35] we obtain 


Pf 


Pi 


0 ^ ^AE-Afi 



Ao 


( 8 ) 


Comparing with the equilibrium expression (8) there is the extra term of the second 
brackets, which becomes 1 for Ap, = 0. The precise effect of this extra term in sensing is 
discussed in the next section. For this discussion it is convenient to dehne the effective 
dissociation constants 


Ko = — and Ki = —e^^. 
a a 


( 9 ) 


3. Integrated sensitivity 

A key observable in sensing is the sensitivity 

r) r)P 

P(c, AE) = -4^:—{a)c = -4 ^ 


9 Inc 


9 Inc’ 


( 10 ) 
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Figure 2. Increase in sensitivity i?(c) through nonequilibrium driving. The full blue 
region under the curve corresponds to the equilibrium contribution ln{Ki/K q), while 
the striped region highlights the additional nonequilibrium enhancement, which is equal 
to the driving affinity Ay. The parameters are set to Ki = 1 /Kq = e^'^, a = e®, and 
Ay = 2.5. 


which is the response of the average activity to small changes in ligand concentration. 
It is convenient to rewrite the sensitivity as 

fl(c,Z\£) = 4P„Pi^ln:^<^ln^, (11) 

a In c Pi a In c Pi 

where the inequality comes from PqPi <1/4. From Eq. (8) it follows that the upper 
bound on the right hand side of Eq. (11) does not depend on AE. Hence, for given c 
there is an optimal conformational free energy difference that maximizes the sensitivity 


AE*{c) = Z\/i +In 


1 + ^, 


+ In 




( 12 ) 


which is obtained from Eq. (8) with PqPi = 1/4. A free energy close to this optimal 
value can be achieved through an adaptation system that uses the methylation levels 
to adjust AE in accordance with the external concentration [13]. From now on we set 
AE = AE*{c) and denote this maximal sensitivity by P(c) = P(c, AE*{c)). 

For the equilibrium case, expressed in Eq. (2), the sensitivity becomes 
, c c 


P, 


eq 


C = 


Ko C + Ki 

Whereas out of equilibrium, with the ratio of probabilities in Eq. (8), we obtain 
P(c) Peq(c) “b Pneq(c), 

with 

Pneq(c) 


(13) 


(14) 


(15) 


c+{^) c+{^) c + Ko c + K,' 

Therefore, the effect of adding a driving affinity Afi to the single receptor is to increase 
the sensitivity by Pneq(c). Particularly, the sensitivity of this single receptor out of 
equilibrium is equal to the sensitivity of a equilibrium model that has a second binding 
site with dissociation constant Ka, as given by (9). This situation is represented in 
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Figure 3. Effect of parameter a on the sensitivity R(c). The parameter a is indicated in 
the figures, Ki = 1/Kq = e^ ®, and Ay = 2.5. For a = e^ the driving affinity increases 
the range for which the sensitivity is non-negligible, while in the other two cases we 
can see a clear increase in the sensitivity in the equilibrium range Kq < c < ifi, with 
a = e^'^® corresponding to the optimal increase in the equilibrium range. 


Fig. 2, where we show the equilibrium contribution to sensitivity peaking between the 
concentration range Kq < c< Ki and the nonequilibrium contribution peaking between 
the range Kq < c < Ki. Calculating the maximum of R^eqic), with c as the optimizing 
parameter, we obtain the inequality 

gZ\/i/2 _ ^ 

-Rneq(c) < ^ ‘ 

We note that an enhancement on sensitivity due to an nonequilibrium driving affinity 
has been shown in [9] (see also [10]). 

As a hrst main result, we obtain that the integrated sensitivity I has the following 
simple relation with the driving affinity Ap, 

/ OO 

d{\n c)R{c) = Ap + In (17) 

-oo Ro 

where we used Eqs. (13), (14) and (15). This result for the integrated sensitivity provides 
a precise quantihcation of the effect of free energy dissipation on sensing. The integral 
represents the area under the curves in Fig. 2 where the equilibrium contribution Req{c) 
yields \n{Ki/K q) and the nonequilibrium contribution rises the area under the curve by 
Ap. 

From the expressions (13) and (15) it follows that Req{c) < 1 and i?neq(c) < 1, 
respectively. The effect of the driving affinity on the sensitivity is twofold: it can 
increase the concentration range for which the sensitivity is non-negligible and it can 
increase the sensitivity in the equilibrium range Kq < c < Ki. The nonequilibrium 
enhancement can even lead to R{c) > 1 within this region. 

The influence of the parameter a as determined by Eq. (15) on these two effects is 
shown in Fig. 3, which indicates that there is an optimal a for which the effect of Ap 
is mostly to increase the sensitivity in the equilibrium range Kq < c < Ki. To quantify 
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^Sa', = r d(\nc)R„^(c) 
JlnKo 


(18) 


We point 

-^neq(c) ^ 


out that due to Eq. (17) and Iko,Ki — due to 

1. Maximizing this integral with respect to a we obtain 

2n 


rneq,opt — rneq _ i 

^Ko,Ki - ^^^Ko,Ki - 



Aii,/2 / Kq 
y Ki 

QAfi/2 _|_ 


(19) 


where the maximum is obtained for a = \JKq/Ki exp(Z\/i/2), which leads to 
dissociation constants that satisfy KqKi = KqKi. Hence, expression (19) provides the 
optimal sensitivity enhancement due to the driving affinity Api within the equilibrium 
concentration range Kq < c < Ki for given Ap, Kq, and Ki. We note that the effect 
of increasing the sensitivity beyond the equilibrium range can represent an important 
advantage for the cell. This increase is quantihed by / — \\i{Ki/K q) — ^ 

more quantitative relation could arise from studying the sensitivity integrated over 
some concentration range of interest. Our choice in Eq. (18) is motivated by the fact 
that convenient for the analogy between nonequilibrium sensing and kinetic 

proofreading. 


4. Analogy with kinetic proofreading 

In this section we establish an explicit analogy between sensing with nonequilibrium 
receptors and kinetic proofreading, with the integrated sensitivity in Eq. (18) playing 
the role of the error reduction due to dissipation in kinetic proofreading. 

4.1. Kinetic proofreading 

The model for kinetic proofreading is illustrated in Fig 4. Two substrates S = R, W, 
with R being the “right” substrate and W the “wrong” substrate, can bind to the enzyme 
E. In equilibrium, the substrate R is copied to a template with higher probability due 
to a free energy difference AF. Specihcally, this free energy difference between state 
EW and ER leads to an equilibrium error 

Ceq = exp(-AF), (20) 

where the error is dehned as the ratio between the probability of writing W and the 
probability of writing R to the template [22,23]. 
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Figure 4. Model for kinetic proofreading. The difference in the transition rates for the 
two cycles is in the free energy term e^^ in the transition from EW to E and from 
EW* to E, which is related to the higher free energy of EW in comparison to ER. The 
rate at which information is written w is assumed to be small compared to the other 
transition rates. 


In the kinetic proofreading scheme phosphorylated forms of the substrates are 
added, leading to the additional states EW* and ER*. The transitions in Fig. 4 involve 
phosphorylation reactions 


ES + ATP ES* + ADP 

(21) 

and dephosphorylation reactions 


ES* E + S + Pi. 

(22) 


If a cycle E —ES —ES* —)■ E is completed, one ATP is consumed and ADP + P, is 
produced, leading to a free energy consumption of Api = Patp — hAOP — hPi- The specihc 
transition rates are shown in Fig. 4, where fc_, /c+, 7 , and 0 are kinetic parameters. 
Moreover, w is the rate at which the substrate S is written to the template, which we 
assume to be much slower than the other transition rates, i.e., we assume the limit 
w —0. 

The error is given by 

^ Pew* ^ f (e~^^ + 0)7 + 0^- \ ( {I + (f))l + 

-Per* + (;/)) 7 e^^ + 0fc_e2"^^/ \ (1 + 0)7 + /’ 

where Pew* and Per* denote the stationary probabilities of states EW* and ER*, 
respectively. Hence, as hrst observed by Hopheld and Ninio [22,23], with energy 
dissipation the error can be smaller than Ceq. The maximal error reduction e/ceq = 
takes place for an appropriate choice of the kinetic parameters and the formal limit 
Ap, —00. 

4.2. Non-equilibrium sensing vs. kinetic proofreading 

The minimal error Copt for hxed free energy difference AF and driving affinity Ap, that 
is obtained by optimizing e in Eq. 23 with respect to the kinetic parameters, is given 













error reduction 


Nonequilihrium sensing and its analogy to kinetic proofreading 


9 



(a) 



(b) 



Figure 5. Kinetic proofreading (left panel) versus nonequilibrium sensing (right panel), 
(a) Error reduction e/ceq and optimal error reduction Copt/ceq, as given by (24), as 
functions of the affinity Ay. The green dashed lines indicate the asymptotically reached 
bounds Copt/eeq > e“^^ and Copt/eeq > . (b) Nonequilibrium contribution to the 

sensitivity integrated in the equilibrium range its optimized value 

as given by (19), as functions of the affinity Ay. The green dashed lines indicate 
the asymptotically reached bounds I’kqKi — I'kqKi — ^rL{Ki/Ko). (c) The 

discriminatory index as a function of the free energy difference AF'. The lower 
panel shows r'opt = —dAP^Tn-^optiAF, Ay), i.e., the discriminatory index associated 
with the minimal error (24). The highlighted areas illustrate relation (27). (d) The 
sensitivity R{c) and the sensitivity i?(c) |opt, which is associated with The 

highlighted areas illustrate relation (18). Parameters are set in the following way: 
AF = \n{Ki/Ko) = 5 with Ki = I/Kq in (a) and (b); Ay = 2.5 in (c) and (d); in (b) 
and (d) I’koKi obtained from (18) with a = e®; in (a) and (c) e/ceq is obtained from 
(23) with fc_ = 7 = 1, ^ = 10“"^. The dotted vertical line in (a) and (b) indicate the 
affinity Ay = 2.5. The dotted vertical line in (c) indicates AF = 5. 
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by [26,27] 


"-^{AF,Ap)=e-^^ 

Ceq 


Afj. AF 

6 2+62 

Afi AF 

6 2+6 2 


(24) 


Sinc6 this function is bound6d by 6 and by 6 th6 following inoquality holds [26], 


e > Copt > 6 xp(-Z\F - Afi). (25) 

Comparing oxprossion (24) for tho maximal orror roduction in kinotic proofroading 
with 6 xpr 6 ssion (19) for tho maximal incroaso in tho intogratod sonsitivity in 
nonoquilibrium sonsing, a quito transparont analogy aris 6 s, as shown in Figs. 5a and 
5b. Both 6 xpr 6 ssions aro tho samo with tho incroaso in sonsitivity in tho oquilibrium 
rang 6 boing analogous to — ln(eopt/eeq) and tho ratio of tho dissociation constants 

Ki/Kq being analogous to 6 "^'^. Whereas in kinetic proofreading a driving affinity Ap 
decreases the error, in nonequilibrium sensing Ap, increases the integrated sensitivity in 
the equilibrium range Kq < c < Ki. 

A recently introduced quantity in kinetic proofreading is the discriminatory index 

| 32 | 

v(AF) = --^ln€, ( 26 ) 

where e is given by (23). Due to egq = the discriminatory index is z/gq = 1 in 

equilibrium, with a larger index > 1 requiring energy dissipation. We can rewrite (26) 
as 

.AF 

-In—=/ dAF'[iy{AF')-l]. (27) 

^eq Jo 

Comparing this equation with (18) we observe that the discriminatory index is analogous 
to the sensitivity R{c), with v{AF') — 1 being the nonequilibrium contribution. In Figs. 
5c and 5d we show the comparison between discriminatory index in proofreading and 
sensitivity in nonequilibrium sensing. Murugan et ah [32] have shown that the integral 
from —cx) to cx) of lo^AF) — 1 can be equal to Ap. This result is equivalent to our 
equality (17). 

5. Effect of the occupancy of the receptor on phosphorylation and 
dephosphorylation rates 

We now generalize the model from Fig. 1 to include phosphorylation and 
dephosphorylation reactions for both 5 = 0 and 6 = 1. With this generalization there are 
two links for the vertical transitions in Fig. 6 a. These reactions happen with transition 
rates 

@ + ATP [I] + ADP ^® + ADP + Pi, 

K*’ Ul*’ 


( 28 ) 
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Figure 6. Four-state system with phosphorylation and dephosphorylation for both 
& = 0 and 6 = 1. (a) Full model with two links for the vertical transitions. The dashed 
links indicate transition rates that are zero in the model of Fig. 1. (b) Total rates as 
given by (32) and (33). The coarse-grained entropy production (43) is calculated with 
these total transition rates. 


where b = 0,1. For thermodynamic consistency, the following constraints must be 
fulhlled: 


for b 


and 


0 , 1 , 


Ap = In 


In 




In 




b b 

(29) 

K_UJ_ 


= z\/i +In 

Ao 

(30) 

A 1 

= A/i-ln—, 

Eo 

(31) 


where we used the free energy (1) for the second and third equations. Whereas the 
presence of two links is important for calculating the rate of dissipation in this model [33] , 
for the purpose of calculating the stationary probabilities we consider the total transition 
rates from inactive to active 

rX = KX + ut (32) 


and from active to inactive 

e^^rt = rff +a;+, (33) 

which are indicated in Fig. 6b. We choose the rates nfk and to be proportional 
to e^^, which leads to independent of AE. Assuming that the binding/unbinding 
transitions are much faster, the ratio of stationary probabilities (8) for this more general 
model becomes 


Pi 



'’v + ) 





C 

Ko 


c 

Ki 


X(c)- 


(34) 
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As in section 3, the sensitivity (11) is maximized for Pq/Pi = 1, which is achieved at 


ZiB'M = 111 I 1^ I -liix(c), 


(35) 


The equilibrium contribution to the maximal sensitivity is given by (13), while the 
nonequilibrium contribution is 


-^neq(c) — 


^ 1 / 3 1 x'(c) 

Inx(c) = -- 


5(lnc) ^ c x(c) 
The integrated sensitivity then becomes 


(36) 


/ = 


d(lnc)[i?neq(c) + Req{c)] = lu ^ + lu 


= In 


^0 ^ 1 


^0 1 


(37) 


For Api > 0, from Eqs. (29), (32), and (33), we obtain 


a;' 

^ < 




QAEj.lp 


(38) 


With these inequalities, we obtain that the integrated sensitivity (37) is bounded by 


J = In 


yt 0 yi 1 


^ 0^1 / 


< In 


K+cai Ki 

= In ^ 

Khcai Kq 


(39) 


where we used (30) in the last equality. As shown in Appendix A, this inequality can 
also be generalized to an arbitrary number of binding sites. 

The influence of how the occupancy of the receptor affects the reaction rates for 
activity on the relation between the integrated sensitivity I and the driving affinity Afi 
can be seen with the following examples. First, if we choose transition rates satisfying 
the relation 


Korf 

K,rl 




(40) 


the function y;(c) in (34) becomes independent of c. In this case, from (36) we obtain 
-Rneq(c) = 0, whicli implies / = ln(iFi/iFo)- Hence, it is possible to have a dissipative 
model with ATP consumption that has the same sensitivity as the equilibrium case. 
Second, we consider the case where phosphorylation happens only if the receptor is 
bound and dephosphorylation occurs only if the receptor is unbound, which is the 
opposite of the model in Fig. 1. In this case n\= uj\ = 0, leading to 


I 1 1 


Z\/i, 


(41) 


where we used (39) and (31). This result shows that the integrated sensitivity can also 
decrease with energy dissipation. The regime for which the integrated sensitivity is 
decreased by Ap, is equivalent to an anti-proofreading regime recently studied in [32]. 
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A more precise analysis of the relation between / and energy dissipation can be 
achieved by considering the entropy prodnction a [33]. For the present model this 
entropy prodnction is the rate of ATP consnmption. Using the stationary probability 
Pa,b we define the probability cnrrent Jb = — Pi,bi^-- With this cnrrent the 

entropy prodnction can be written as 


O' — {Jq + (42) 

by nsing the fact that a is a snm of cnrrents mnltiplying cycle affinities [33] . The energy 
dissipation a is non-zero whenever Ap, ^ 0. Besides a, we can consider a coarse-grained 
entropy prodnction d, which does not take into acconnt the two channels for the vertical 
transitions in Fig. 6a: it is obtained by considering the single links with rates r\ in Fig. 
6b. This coarse-grained entropy prodnction provides a lower bonnd on the fnll entropy 
prodnction, i.e., a > a [36]. For the model in Fig. 6b, 


— .^0 


a = JA. 



Po,o-e^MPi.o 

and 


' 0 1 \ 


M = In 1 

rfrt \ 

— In 

O i 


Ko 

is an effective affinity. From relation (39) we obtain 


(43) 


(44) 


A = I-\n^, (45) 

Ao 

which is the noneqnilibrinm contribntion to the integrated sensitivity I. The effective 
affinity associated with the coarse-grained entropy prodnction determines three different 
regimes for noneqnilibrinm sensing. For ^ > 0 the integrated sensitivity is increased in 
relation to its eqnilibrinm valne. If ^ = 0, which implies d = 0, the energy dissipation 
has no effect on sensitivity. If ^ < 0 then the ineqnality d > 0 implies J < 0. In this 
last regime energy dissipation decreases the integrated sensitivity. 


6. Conclusion 

We have characterized the enhancement of sensitivity by a noneqnilibrinm driving 
affinity that arises from ATP hydrolysis in the chemical reactions involving an activity 
change. For the single receptor model from Sec. 2, the integrated sensitivity I was 
shown to have a simple relation with the driving affinity in Eq. (17). We have shown 
that a dissipative sensing model can lead to both an increase in the concentration 
range for which the sensitivity is non-negligible and an increase in the sensitivity in the 
eqnilibrinm range. The second effect is qnantified by which is defined in Eq. 

(18). 

We have shown that noneqnilibrinm sensing is eqnivalent to kinetic proofreading, 
with the analogons parameters, observables and relations snmmarized in Tab. 1. Most 
prominently, while in noneqnilibrinm sensing a driving affinity leads to an increase in 
the sensitivity integrated over the eqnilibrinm range, in kinetic proofreading a driving 
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Nonequilibrium sensing 

rneq 

R 

ln(iFi/iFo) 

I < ln{K,/Ko) 

-|- Ap, 

Kinetic proofreading 

-ln(e/eeq) 

V 

AF 

e > exp(-Z\F 

- 


Table 1. Nonequilibrium sensing compared to kinetic proofreading. 


affinity decreases the error. In kinetic proofreading the equivalent of sensitivity is the 
discriminatory index introduced in [32]. 

The influence of the occupancy of the receptor on the phosphorylation and 
dephosphorylation rates is of fundamental importance for the relation between 
integrated sensitivity and the affinity. As we have shown in section 5, it is even 
possible to have a regime where energy dissipation leads to a decrease on the integrated 
sensitivity, which is analogous to the anti-proofreading regime from [32]. 

Our results demonstrate that measurements of the integrated sensitivity 
could unveil how the occupancy of the receptor affects the phosphorylation and 
dephosphorylation rates. It is certainly intriguing to speculate whether real chemotaxis 
networks evolved in such a way that this influence optimizes the enhancement of 
sensitivity due to energy consumption. 

Acknowledgments 

D. H. acknowledges helpful discussions with J. M. Sanchez. 


Appendix A. Generalization to an arbitrary number of binding sites 


In this Appendix, we consider a generalization of the model studied in the main text 
for an arbitrary number of binding sites. In this case, the variable b takes the values 
b = 0,1,..., N, where N the number of binding sites. The transition rates for this more 
general model are shown in Fig. Al. Rates involving a change in activity are given by 
(32) and (33). Rates related to a change in the occupancy of receptor must fulhll the 
generalized detailed balance relation with respect to the free energy (1). We set these 
rates as follows. The binding rate from 6 to 6 -|- 1 is Wbb+i = 7b(A^ — b)c/Ka, where the 
factor N — b comes from the fact that there are N — b free receptors for the ligand to 
bind; the unbinding rate from 6 to 6 — 1 is ~ Ibb, where the factor b is related to 

the number of bound receptors. 

We assume that the binding events are much faster than changes in activity. In 
this case, the stationary conditional probability reads 


leading to 


P{b\a) 





Pi 





N 

X(c), 


(A.l) 


(A.2) 
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Figure Al. Generalization of the single receptor model in the main text to an arbitrary 
number of binding sites N. The rates r^j_ and are defined in (32) and (33), 

respectively. 


where 


y-Af /AfW_^^6 6 

j-- 3j 

This expression generalizes (34) to the case of N binding sites. Following the same 
procedure from section 5, similarly to (35) the sensitivity (11) is maximized for 


AE*{c) = iVln 


i±iL 

1 + t 


- Inx(c), 


where, similarly to (36), the nonequilibrium contribution to sensitivity becomes 

^ ^ Ix'(c) 

nn —l™x(c) =-yw. 

a(lnc) c x(c) 

Hence, the integrated sensitivity reads 


-^neq(cj 


I = 


d(ln C) [i?eq(c) + i?„eq(c)] = iV lu ^ + lu 

Kq X( 0 ) 




(A.4) 


(A.5) 


(A. 6 ) 


The transition rates for this model must fulhll the constraint 

f nP. ujP \ iFi 

= N\Yi — + Ap. 

\tx_uj_: J Kq 

From the inequalities (38) we finally obtain 

I < A^ln ^ + Afi, 

Aq 

which generalizes inequality (39) to the case of N binding sites. 


(A.7) 


(A. 8 ) 
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